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Abstract 

We show that correlation functions of the Z„-Baxter model in the principal regime satisfy 
a system of difference equations. We obtain the spontaneous polarization of the Z^-Baxter 
model as a solution of the simplest difference equation. 
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1 Introduction 



It is widely recognized that systems of holonomic difference equations commonly appear in 
integrable quantum field theories and solvable lattice statistical models. Smirnov found 
that form factors in two dimensional integrable massive models satisfy difference equations. 
Frenkel and Reshetikhin showed that correlation functions of the intertwining operators 
of the quantum affine algebra obey difference equations, and that the resulting connection 
matrices give elliptic solutions P, 3 to the Yang-Baxter equation (YBE). In |^ it was 
proposed that correlations of the XXZ model in the anti-ferromagnetic regime can be for- 
mulated in terms of the trace of products of the intertwining operators of Ug{sl2)- The higher 
spin analog of the XXZ model [0] and the higher rank generalization Q were also studied. 
The integral formula of correlation functions of the XXZ model is given in 0. That of the 
matrix element of product of intertwining operators of t/g(sl2) for arbitrary level is obtained 
in 
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It was remarked in that even if the trigonometric /^-matrices is replaced by the elliptic 
ones, integrability of difference equations survives. Jimbo, Miwa and Nakayashiki [|11|] actually 
obtained difference equations for correlations of the eight- vertex model in the principal regime. 



In the present paper we generalize their results to the case of the Z„-Baxter model |]T2 . 

The Z„-Baxter model is a vertex model on a two-dimensional square lattice C such that 
the state variables take on Z„-spin. Each oriented line of C carries a spectral parameter 
varying from line to line. We assign Z„-valued local states on edges. Set 















I 





be a local Boltzmann weight for a single vertex with bond states i,j, k, I G Z^. Arrows denotes 



orientations of lines. In this case the weights are given as follows ||T3 

h{z)e^'-''\z + w) . , 

Riz)^ = \ ^0--^)(.)^(--^-)(^) ^i^ + k = J + l,^odn, ^^^^^ 
otherwise, 



where 



:= y exp L^/^(m +- - ^-fnr + 27rv^(m + - - + - 
I 2 n 2 n 2' 
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n— 1 n— 1 



hiz) := n ^^^H^)/ n ^^^Ho), 

and w 7^ mod Z + Zr, is a constant. Now we assume that 0<t<q<u<l, where 
t := exp(7r-\/— Ir), g := exp(7r-\/^^w), and u := exp{—TT\^—lz). Then the elhptic theta 
functions are expressed in terms of the product series 

= y^c^i/2tn(l/2-i/„)2^-l+2,7n(t2„. ^2n)^(^2j^2. ^2n)^ (^2{n-i)^-2. ^2„)^^ 

where 

oo oo 

(a;gi,---,gfc)oo:= I] •■■ 11 (l-aqr-'-qD. 

mi=0 mj.=0 

Following Baxter we call such domain of parameters the principal regime. Note that ( p,.l| ) 



is weights of the eight- vertex model when n = 2. 

The plan of this paper is as follows. In section 2 we summarize our result on the unitarity 



and the crossing symmetry of the Z„-Baxter model [[151 , p!q| . In section 3 we introduce the 
normalized i?-matrices by the partition function per site in the thermodynamic limit. In 
section 4 we show that correlation functions defined by using this normalized i?-matrices 
satisfy g-difference equations. In section 5 solving the simplest one we obtain the spontaneous 
polarization. By setting n = 2, our result reproduces that of [0. We can also get Koyama's 



result p by taking the limit t ^ 0. In section 6 we give some remarks. 



2 Unitarity and Crossing Symmetry 

Let V = and {vi}ieZn be the standard orthonormal basis of V. Then the i?-matrix R{z) 
whose (i/c, j7)-element is R{zyjl gives a linear map V ® V ^ V V 

i,k£Zri 



The Z„-Baxter model has the Z„-symmetry 

(i) R{zyjl = 0, unless i + k = j + I, mod n, 

(ii) i?(2;)}tp,?+p = Riz)fi, for every i, j, k, I and p E Z„. 

In terms of two linear map in V 



{2.1] 



gvi = u'vi, hvi = Vi-i, (2.2) 
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where, uj = exp(27rv— T/Ti), (|2.1|) can be rephrased as follows: 

R{z){g®g) = {g(^g)R{z), 



(2.3) 

R{z){h®h) = {h®h)R{z). 



Owing to (|2.3| ) the Z„-Baxter model is called the g[„-invariant model. Furthermore, the 



i?-matrix satisfies the Yang-Baxter equation 

Rl-a{Zl - Z2)Rl-'i{Zi - Z-i)R2-3{z2 - Zs) = ^2;3(2;2 " 2;3)^1;3(^1 " 2;3)i?i;2 (^1 - Z2) , (2.4) 

where Ri-j{z) denotes the matrix on V'^^, which acts as R{z) on the i-th and j-th components 
and as identity on the other one. In what follows subscripts of any linear operator signify the 
components on which it acts. 

Let us reformulate this lattice model in terms of the scattering theory of elementary 
particles. Consider the system consisting of n kinds of particles with equal mass and let i* 
stand for the antiparticle of the particle i. In a two dimensional integrable field theory, both 
rapidities and Z„-current preserve before and after any two-body collision. Let V = and 
Vz be a copy of V with a rapidity (spectral) parameter z. Then R{zi — Z2) can be regarded 
as a linear map i?^^i'^^2 acting on V^^ V^a- is due to the Lorentz invariance that i?^^i'^^2 
depends only upon the difference of rapidities zi — Z2. We also set R^^i^^^i = Pi?^^i'^^2^ where 
P is the permutation V^^ (g) V^^ V^^ ® V^^ . In this notation, YBE (|2.4|) reads as follows: 



(/ ® i?^-i'^-2)(i^v'zi,v.3 ^ /)(/ ^ R^^^y^s) = (^vi^y.g /)(/ Rv^iy^3)[Ryny^2 ^ /), (2.5) 

ClS Sb linear map Vz-^ ® Vz2 ® Vz^ Vz^ Vz^ ® • 

Next we extend this formulation for arbitrary spaces K and L of tensor product of some 
V^. For K = Vz^® ■ ■ ■ ®Vz^. and L = \4' ® ■ • • C?) V^/, it is very natural to define a linear map 



i?^'^ ■.K^L^L^Ka.s follows P, 115 



t>Ky, f^z^yz' fy^kyz' 

n .— ni.2 ■ ■ ■ J^k;k+i ■> 

tdK,L ._ id I D ^1 

YBE holds for R^'^ by virtue of YBE for R^'^ (U) 

(/ ® i?^'^)(i?^'*' ® /)(/ ® i?^'"^) = {R""^' ® /)(/ ® i?^'^)(i?^'^ ® /), (2.6) 

linear map K®L®M^M®L®K. 
Now we summarize the unitarity and the crossing symmetry of the Z„-Baxter model. 
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'1) The unitarity The unitarity or the first inversion relation |T^, 16 1 is given by 



h{zi - Z2+ w)h{-zi + Z2 + w] 



(2.7) 



(2) The crossing symmetry Let V* be the dual space of V and {v*}i(=.z„ be the dual 
basis of {vi\i(^Zn- Then we have the isomorphism C -.V* ^ A'"^^(V^) 



cv:= E 



^\■.■■■,^n-l 



n- 1)! 



(2.8) 



where e*^ is the n-th order completely antisymmetric tensor. 

The i?-matrices corresponding to the collision between a particle and an antiparticle is 
given as follows: 



= (C®/)-l/2^-l'^-2 + (n-l)»®-^^-2+-(/® C), 
= (/OC')-l/?^n+(-l)-®-®^-l+-"^-2(C(g)/). 



The matrix elements of these are 



(2.9) 



)ji* 1 



i* ,k 



which meet the crossing symmetry [TB] 



p=2 



h{w) 



]-[ K-zi + Z2- pw) 



h{w) 



From (|2.71 )and (|2.10|) , we have the following second inversion relation 16 



;2.io) 



The i?-matrix corresponding to antiparticle-antiparticle scattering is also defined by 



(2.11) 
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3 Partition Function and /S-Matrices 



Let k{z) be the partition function per site in the thermodynamic hmit. With the help of two 
inversion relations ( p.7|) , (|2.11|) and Baxter's corner transfer matrix method fl^, we can get 
the functional equations for k{z): 

h{z + w)h{~z + w) 



K[Z)K[—Z] 



h\w) 

h(—z)h(z + nw) 
K[z)K[—z — nw) = - 



k{u) 



(^2+2„^-2. q^n^^l^u^- ^ g2n)^(^2^-2+2n^-2. ^2^ ^2n)^(^2^2. ^2^ g2n)^ 

For later convenience, we define S^'^ = PS^'^ by 

n-l 

= \{k{zi- Z2-pw)-^R^^^^^*^, 

p=l 

n-l 



p=l 

n-l 
p,q=l 



Then the unitarity and the crossing symmetry for S ' hold: 



The first one is obvious. The last two follow from 

n— 1 n—1 



K[ —Z 



-1 



n '^(^ - = n 



h{—z + jvj) 
h{w) 



(3.1; 



Hereafter k,{z) is often denoted by k{u) through the relation u = exp(— tt^— I2;). 
In the principal regime using ( |1.2D the following expression solves ( |3.1| ) [|T^ 

K{u) = (3.2) 

where 

iq'u'; t\ g2«)^(g2n^-2; t\ q'-)ooit'q-'u'; t\ q^-)^{fq'-u-^- 1\ q^^U 



(3.3) 



(^^^i'^^2)p = (^^^2.v.i)g, (3.4) 
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which can be checked by exphcit calculation. The relations (|3.4| ) are visualized as 




I* 



Z2 



I* 





'Z2 




Zl 


I 










.Zl 


I 


Z2 



where z = z + nw. Note that z (z*) in the above figure represents the space (V*), while i 
(i*) denotes Vi (v*). 

Remark Speaking in terms of Young tableau, y = □, the fundamental representation 
of sin, while V* corresponds to n — 1 vertical D's, the space of antisymmetric tensors in 
Thus both V and V* are indispensable in our formulation. For n = 2 |jTl| from (|2.8|) and 
]9D, one can do without V* by identifying v* = (— V* = Vz+w 



4 Difference Equations for Correlation Functions 

In the principal regime the Boltzmann weights of the types S'(2;)*^J dominates the others. 
Thus in the low temperature limit t,q ^ 0, only the configuration such that the spin variables 
take the same value in the direction from NE to SW and increase by one in the direction from 
NW to SE, is possible. We call it a configuration of the ground states. There are n ones 
labeled by m G Z„. In what follows, we fix one of them (say, m) and define all the correlation 
functions in terms of the low-temperature series expansion (i.e. the formal power series of 
t and q). Then the lowest order of them comes from the m-th ground state configuration. 
Furthermore, any finte order contribution is derived from the configurations which differ from 
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that of the m-th ground state by altering a finite number of spins. It is equivalent to taking 
the GNS representation obtained from the m-th ground state (m-th GNS representation) as 
the Hilbert space. It is expected that the correlation function defined in such a way is an 
analytic function which has a finite convergence radious if there exists the phase transition 
at a finite temperature. 

Let us consider the probability 

Zi Zn 



(\ h,---,iN 
Ai,-/^ config 



in the inhomogeneous Z„-Baxter model. Here m specify the boundary condition such that the 
spin on the reference edge equals to m in the ground state configuration. The reference edge 
is the next left to the one with the spectral parameter Zi. The symbol of the sum denotes the 
statistical sum over the m-th GNS representation. If we assign the weight S to each vertex, 
the statistical sum gives just the probability in the thermodynamic limit. 

The YBE, the unitarity and the crossing symmetry permit the following manipulations 



In 
In 
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(m+ 1) 







{m - 1) 



where z — z + nw. 
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Thus we conclude that the probabihty P^"^^ 




is given by the following 



correlation on the dislocated lattice 



config t ft 



If we set z,, = zl,i„ = < v < N) in the above expression, we can get the probability 
such that the spins on N successive vertical spins located in the same row take the values 
ii, ■■■ ,iN- 

We also define the correlations 



config 



Ci C2N 

where the set of Cj's is a permutation of z[, - ■ ■ ,z'jq,z{, - ■ ■ ,zlj, and the set of ij's is that of 
i'l-i ■ ■ ■ ) ^'n-: ill' ' ' 1 i*N- Consider the V(^^® ■ ■ ■ ® V(^^j^-Ma\Med correlators 

FM(Ci,...,C2iv) = E v,,®---®v,,,F^^){Cir--.C2N\„-,.,, (4.1) 
where Vz* — V* and Vi* — v*. 
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The correlators ( |4.1|) have the S'-matrix symmetry: 

F^^'H- ■ ■ , 0+1, 0, ■ ■ ■) = SJ^^'^^F^\. . . , 0, 0+1, ■ ■ ■)• (4.2) 
Furthermore, it satisfies 

F^"'^'\C2, ■ ■ ■ , C2N, Cl) = P2N^l;2N " " " Pl;2F^"'\Cl, C2, " " " , C2n). (4.3) 

In the left hand side we should take + from the double signs if Ci = z^, (1 < 3z/ < N), — if 
(i = 2;*, (1 < 3/i < A^). The symbol ( means z + nw if ( = z, {z + nw)* if C = By setting 
n = 2 eqs. (^^) and ( [4.3|) give those of the eight- vertex model In these are two of 
axioms that the form factors of integrable massive models should obey. As for the derivation 
of them, see 

From ( [4. 21 ) and ( [4. 3D we obtain the difference equation 

F^-^'K- ■■,(„■■■) = (5^5^^"^)-^■■■(5^|f^")-^x 

x^:^'^^...s;;!^'^^FM(...,o,---), 

where the interpretation of the left hand side is the same as that of ([4.3|) . 



5 Spontaneous Polarization 



In this section let us concentrate to F^"^\z', z*). Since it depends only upon u = exp(—n^/^(z— 
z')), we denote F("^)(z',z*) by F^'^^u). Set 

n-1 



G('")(w):=^^"^'F«(n)oo.. (5.1) 



1=0 

Then we have 

, ' = Vu;™(^--^)(^^^"^^*)°°:, (5.2) 

where we use the Z„-symmetry F^™'\u)jj* = F^^'^^p^u) j+p(^j+py . Let us solve ( |5.2| ) under the 
condition lim^^g^o ^^"^(1) = which follows from that F^'^^z' , z*)jj*\z=z' SJ''^^ in the 
low temperature limit. 

Using the crossing symmetry for S'-matrix, we get 

^{n-2)/n(^2n.^2n)2^ ^^2 . ^2)^ ^^2^-2 . ^2)^ + 



^ Riu) (t2;t2)2^ („2.^2n)^(^2n^-2.^2n)^ QU) ( 
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It follows from the transformation properties for the theta functions, 

^^^_^^ e^^\z-z' + w) ^ _^ K{z - m)/n + w) Yli^^ h{{-z + l)/n) 

{q^;t^)oo{t'q-^;t^)oo 

where we use ( |1.2|) at the second equality. Therefore (|5.2|) reduces to 



Note that the relation 



(5.3) 



(p{uq- 



where 

^(u) := giuq-/')g{u-'q-/'), g{u) : 
Then we obtain 



(g^^M-^; t^, g2n^ ^2n)^^^2g3n^-2. ^2^ ^2n^ ^2n)^ 



(g2+n„2. ^2^ ^2n^ ^2n)^ (^2g-2+n„2. ^2^ g2n^ ^: 



2n^ ■ 

oo 



^ ^ <^(1) (t2;t2)2^ (g2^™n2/";g2)^(^2^-m„-2/n.^2)^- ^ ' ) 

The uniqueness of this solution is ensured under the assumption of the analyticity. 

If we define E^^'^u) = E]Zq F^""Ku)jj,, then E'-"'\u) = uj"'G'--^\u). Since it is the 
expectation value of g defined in (|2.2|), the polarization of this model is given as follows: 



(tVt^)oo(t^^-^t^)oo (5.5) 



,2) 

loo 



(t2;t2)2^ (^2^. ^2)^(^2^-1. ^2 

It reproduces the polarization of the eight-vertex model conjectured by Baxter and Kelland 
p!8| and confirmed in |]ll| when n = 2. Taking the limit t — 0, it gives the result of in 
which the one point function of the s[„-analog of the XXZ model is calculated. We can also 
obtain the expectation value of g^ 

(t2;t2)2^(g2^^g2)^(g2^-/;g2)^- ^ ' ) 
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6 Concluding Remarks 



In this article we study the symmetry of the correlation function of the inhomogeneous Z„- 
Baxter model on the dislocated lattice, which coincide the probability of the appropriate 
spin configurations if one tune the value of spectral parameters and spin variables. Such 
correlations can be characterized as solutions of the system of difference equations. By solving 
the simplest one, we obtain the spontaneous polarization. This result includes various ones 



already obtained [0, |T2, P as special cases 



The next problem is to calculate the A^-point function of the Z„-Baxter model. In the 
trigonometric case the free field representation of the g-deformed vertex operator of the quan- 
tum affine algebra and the resultant integral formulae of A^-point functions were obtained 



IU| . Since the elliptic generalization of the quantum affine algebra has not yet been found. 



it is difficult to apply the strategy in ^ |T0| to the present case. Thus it is promising to derive 
and solve the recursion relation of the correlations. 

The critical behavior of this model is also interesting. The six- vertex model with = 1, 
which is related to minimal models of conformal field theories, is the critical limit of the 
eight- vertex model Refs. [H, 01 treat the six- vertex model and its higher spin analog 
in the antiferromagnetic regime (—1 < g < 0) because they are based on the corner transfer 



matrix method ||14[. It is valuable to discuss the above trigonometric models with |g| = 1 as 
the critical limit of the Z„-Baxter model. 
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